Modal dynamics in multimode fibers 
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The dynamics of modes and their states of polarizations in multimode fibers as a function of time, space, and 
wavelength are experimentally and theoretically investigated. The results reveal that the states of polarizations 
are displaced in Poincare sphere representation when varying the angular orientations of the polarization at 
the incident light. Such displacements, which complicates the interpretation of the results, are overcome by 
resorting to modified Poincare spheres representation. With such modification it should be possible to predict 
the output modes and their state of polarization when the input mode and state of polarization are known. 
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1. Introduction 

Multimode fibers have high coupling efficiencies, can op- 
erate over a wide range of wavelengths, and have rela- 
tively low susceptibility to degrading nonlinear effects. 
As a result they are widely used in mid and short com- 
munication ranges, high power lasers and amplifiers, and 
in transporting high optical powers from one location to 
another [1,2]. Yet, while in single mode fibers the po- 
larization dispersion is low and can be significantly re- 
duced either by using polarization maintaining fibers or 
by operating at the principal state of polarization [3-9] , 
in multimode fibers the mode dispersion is high result- 
ing in significant reduction of performance [10,11]. Con- 
sequently, it is important to continuously analyze and 
monitor the dynamics of the modes propagating in mul- 
timode fibers. 

When analyzing the dynamics of modes in multimode 
fibers, it is common to separate the dynamics of the 
transverse modes and the dynamics of the states of po- 
larization [12-15]. Unfortunately, since there is coupling 
between transverse modes with different states of po- 
larization, such a separation is not natural and results 
in coupled equations that are difficult to analyze. Here 
we present a new approach where the complicated dy- 
namics of the modes in multimode fibers are separated 
and represented by two uncoupled equations of motion 
that lead to simple geometric plots on modified and un- 
coupled Poincare spheres. Such plots fully describe the 
modal dynamics in multimode fibers, and can be con- 
veniently used for monitoring, analyzing and predicting 
them. 

We begin by describing our experimental configura- 
tion and present experimental results of the state of po- 
larization (SOP) at the output from multimode fiber as 
a function of wavelength for different angular orienta- 
tions of light with linear polarization at the input. The 
results reveal that the dynamics of the modes and their 
SOP at the output from multimode fibers are relatively 
complicated. Then we developed an analytical model for 
calculating these dynamics, and show how their repre- 



sentation can be simplified to obtain a more convenient 
and intuitive understanding. Finally, based on the results 
from the model, we show how to obtain such simplified 
representation experimentally. 

2. Experimental configuration and results 

The experimental configuration for measuring the out- 
put polarization from a multimode fiber is presented in 
Fig. 1. A linearly polarized light from a tunable fiber 
laser propagates through a half wave plate and a quar- 
ter wave plate in order to obtain light of any desired SOP 
and wavelength at the input. The light is then launched 
into a multimode fiber, and the SOP at the output as a 
function of space, time and wavelength is measured using 
a real-time space-variant polarization measurement sys- 
tem [16]. In our experiments the multimode fiber length 
was 10m, core diameter 20fim and numerical aperture 
0.07. It mainly supported the T-Eqi and TMqi modes 
which in the weakly guided approximation can be ex- 
pressed as orthogonal sets of linearly polarized modes. 
The electric field intensity distributions and polarization 
directions for these modes are presented in Fig. 2. As 
shown, there are four distinct modes of LPn, each with 
different intensity distribution or different orientations of 
the polarization direction. The white arrows at the cen- 
ter of the individual intensity distributions denote the 
direction of the polarizations in each mode. 
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Fig. 1. Experimental configuration for measuring the 
output polarization from a multimode fiber as a func- 
tion of the input wavelength for different input states of 
polarization. 

We started our experiments by verifying that the cou- 
pling into the fiber does not change the mode or the SOP 
of the input. This was done by determining that the in- 
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Fig. 2. Intensity distributions together with the orien- 
tations of the polarization direction for the orthogonal 
sets of linearly polarized modes LPn. The white arrows 
in the center of the individual intensity distributions de- 
note the polarization directions. The black arrows and 
the terms U2 and denote the coupling between the 
different modes. 



tensity distribution of the mode and the SOP at the 
output are the same as those at the input after propa- 
gating through a short 1cm fiber. Then, we launched the 
same linearly polarized light with a uniform intensity dis- 
tribution into the 10m long fiber, and measured the SOP 
at the output while varying the input wavelength from 
1062nTO to 1066nm. We repeated these measurements, 
each time after rotating the linearly polarized light at 
the input by 10°. 

The results are presented in Fig. 3 as a Poincare sphere 
representation of the SOP at the output for different an- 
gular orientations of the linear polarization at the input. 
As evident, the circle of smallest size is obtain when the 
linearly polarized input beam is oriented at 30°, indi- 
cating where the SOP of the output is nearly invariant 
to the input wavelengths. Namely, at an angular orienta- 
tion of about 30° we have a principal mode (PM). As the 
deviation from 30° increases, the circles become larger 
indicating that more spectral dispersion occurs. Unlike 
increasing concentric circles around the principal state 
of polarization location in single mode fiber, here we ob- 
tained increasing circles whose centers continuously de- 
viate from the PM location as the difference between 
the input and the PM increases. We also present the 
projection of the output polarization on the 5*1 axes as a 
function of the input wavelength for two input polariza- 
tion in Fig. 4. As evident, when the input polarization is 
close to the PM the output polarization is less sensitive 
to wavelength variations. 

3. Model and calculated results 

We developed a model for calculating the dynamics 
of the SOP at the output of multimode fibers, in or- 
der to gain a better physical intuition about the evo- 
lution of modes in a multimode fiber and to under- 
stand what causes the deviations from the PM loca- 
tion. We define E.^ and Eout as the input and output 




Fig. 3. Poincare sphere representation of the SOP at the 
output as a function of the wavelength for different an- 
gular orientations of the linear polarization at the input, 
(a) Experimental results; (b) calculated results. The in- 
put wavelength ranged from 1062nrn to 1066nTO for each 
angular orientation. 
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Fig. 4. Projection of the output polarization on the Si 
axes as a function of the input wavelength for two ori- 
entations of the input polarization; asterisks (blue) - 30° 
orientation, close to the PM; circles (red) - 90° orienta- 
tion, far from the PM. 



modes composed from a superposition of the four fields 
Ein.out = [El, E2, E^, Ei]^^^^^^ shown in Fig. 2. The re- 
lation between them is 



(1) 



where j3 is real and U is a unitary matrix. To explic- 
itly define U, we assume that birefringence and twists 
in the fiber are the main sources for coupling between 
the modes. Birefringence in the fiber leads to coupling 
of modes with the same intensity distributions but with 
orthogonal SOP. Accordingly, the coupling between the 
fields El and E2 can be described by a 2 by 2 matrix, as 



U(2)(c^) = 



(2) 



Ui{lo) U2{lo) 
-U*2{uj) U\{uj) 

where ui denote changes in the phase of each mode while 
propagating in the fiber and U2 the coupling strength be- 



2 



tween the modes. This matrix is the same as that for cou- 
pHng orthogonal polarizations in single mode fibers [5]. 
Similarly, since and E4 also have the same intensity 
distributions and assuming that birefringence is essen- 
tially constant across a relatively small fiber core [f7, 18], 
we can describe the coupling between £^3 and E4 with 
the same unitary matrix {/(^^(cj). Twists in the fiber ro- 
tate the modes and leads to coupling between £"1 and E3 
and between E2 and E4 with coupling strength denoted 
as Us. We also assume that coupling between Ei and 
E4 and between E2 and E3 are negligible. Taking bire- 
fringence and twists into account, the overall coupling 
matrix between all four fields is then U(w) = exp(H) 
where H is skew-Hermitian matrix as 



where 



H(a;) = 



■ui(a;) U2{uj) 

















ulioj) 



(3) 



Following the derivation of Poole and Wagner [5] and 
extending it to our multimode fibers [14, 15], leads to 
four principal modes (PMs) where each is insensitive to 
variations in the input wavelength up to the first order. 
The four PMs, denoted as e„ where n = 1, 2, 3 and 4, are 
the solutions of 



[U' - ifc„U] e„ = [H - iknl] Ue„ = 0, 



(4) 



where the prime denotes differentiation with respect to 
frequency, and the four fc„ are 



kn = ±n/ |Wl 



\U2 



|2 + |M3|2±2|U3|V|M1|2 + |W2P. 



(5) 

These four PMs will propagate in the multimode fiber 
with minimal dispersion. 

To analyze the dynamics of modes other than PMs 
as a function of the input wavelength, we differentiated 
Eq. (1) with respect to frequency, and used the definition 
of U to obtain 



(6) 



Using Pauli matrices, a^, Uy and a^, we can decompose 
Hto 

H = {ibia^ + ia2(Jy + ih2(Jx) + {ibso^x + ia^o-y) , 

(7) 

where ai denote the real part of Ui and hi the imaginary 
part of Ui. Equation (7) defines rotation in four dimen- 
sions, one part for left rotation and the other for right 
rotation [20,21]. Now, in order to find analytic solution 
to Eq. (6) we redefine the vectors Ein and Egut to a 2 
by 2 matrices, 



E„ 



in, out — 



El E2 
E3 E4 



(8) 



in, out 



and substituting Eqs. (7) and (8) into Eq. (4) leads to 
the solution of Eout as a function of frequency as 



Eout{t-j) = r2i(a;)£i„ll_R(a;), 



(9) 



cos(|Tb|w) |^sin(|Tfc|w) 
-|^sin(|T6|a;) cos(|r6|a;) 



and 



cos(|ra|w) + 1^ sin(|Ta|w) 

\ra\ 



^sindrjoj) 



COS! 



\ra\ 



sm(|T, 

_ "1 
\r.\ 



^Sm{\Ta\co) 



with Ta = {bi,b2, 02) and n = (0, 03, 63). 

Equation (9) provides an analytic solution of the out- 
put mode distribution and its SOP as a function of the 
wavelength when the input SOP is known. It involves a 
multiplication of two matrices flu and Hl. The matrix 
flfi is composed from ui and U2 so it represents the cou- 
pling between each pair of modes with the same intensity 
distributions but with different SOP whereas the matrix 
Hl is composed only from M3 so it represents the cou- 
pling between each pair of modes with different intensity 
distributions and different SOP. Using Eq. (9), we calcu- 
lated the evolution of the output SOP as a function of the 
input wavelength for different angular orientations of a 
linearly polarized light at the input. The results are pre- 
sented in Fig. 3, with the coupling parameters ui, U2 and 
U3, as fitting parameters. As evident, there is a qualita- 
tive agreement between the experimental and calculated 
results, indicating that our assumptions in developing 
the model are justified. We attribute the difference be- 
tween the calculated results and the measured results to 
small deviations in the polarization measurement and to 
inaccurate coupling parameters. 

4. Modified Poincare spheres representation 

In order to simplify the Poincare sphere representation 
of the SOP dynamics at the output, we resorted to 
modified Stokes parameters Sa = {Si^\ S^^) and 

5, = {Si'\si^\si'^) where 



Si'^ = \Ei\^-\E2\^ + \Es\^-\E4^ 

Si^^ =2\Ei\\E2\cOs{ipi2) + 2\E3\\E4\cOs{>fi34) (10) 

Si^^ = 2|£;i||£;2| sin(^i2) + 2|£;3||E4| sin ((^34) 



and 



5W = |£,|2 



|i?2p-|i?3p-|i?4p 



;.(2) 

:(3) 



2\Ei\\E3\ cos (ipia) + 2\E2\\E4\ cos ((^24) (H) 



S^"' = 2|£i||£;3| sin((^i3) + 2\E2\\E4\ sin (^24) 

where (pij denote the phase difference between E^ and 
Ej. Now, Eq. (6) can be separated into two uncoupled 
torque equations, as 



Sa — Ta Sa 
Sb = Tb X Sb- 



(12) 



Note, the torque vectors depend only on the parameters 
of the fiber and not on the input SOP. 
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Using the two sets of modified Stokes parameters we 
recalculated the evolution of the output SOP as a func- 
tion of the input wavelength for different angular ori- 
entations of a linearly polarized light at the input. The 
results are presented in Fig. 6. As evident, the use of 
the modified Stokes parameters and the separation of 
Eq. (6) into two uncoupled equations, resulted in two 
modified Poincare spheres where the SOPs in each sphere 
form concentric, rather than displaced, circles around the 
PM when varying the input wavelength. The modified 
Poincare spheres representations involve natural modes 
of a fiber. Specifically, Sa represents the modes TEqi and 
TMqi, and Sb represents the modes HE21 and 
These results indicate that when evaluating the propaga- 
tion of light in multimode fiber it is preferable and more 
convenient to resort to our uncoupled modified Poincare 
spheres representation rather than the usual Poincare 
spheres representation [12,13]. 




Fig. 5. Calculated modified Poincare sphere representa- 
tion of the SOP at the output as a function of the wave- 
length for different angular orientations of linearly po- 
larized light at the input, (a) Using the modified Stokes 
parameters Sa ; (b) using the modified Stokes parameters 



Based on the calculated results from our model we 
reevaluated our experimental results in order to obtain 
more convenient and more intuitive Poincare sphere rep- 
resentations of the SOP dynamics at the output from 
a multimode fiber. Specifically, we first multiplied the 
measured light distribution with those of the transverse 
modes shown in Fig. 2. Then we determine the SOP 
at each point across the beam which leads to determin- 
ing the contribution of the fields Ei, E2, E^, and E4. 
Next, we used Eq. (10) to determine the modified Stokes 
parameters Sa to obtain the desired modified Poincare 
sphere representation. For simplicity we present only the 
average SOP (center of circles) for each angular orienta- 
tion of the linearly polarized input beam. The results are 
presented in Fig. 6. Figure 6(a) shows the average SOP 
for different input polarization orientations in the origi- 
nal Poincare sphere representation, same as Fig 3. Here, 
the average SOP is displaced as we vary the input polar- 
ization orientation. Figure 6(b) shows the corresponding 
results on the modified Poincare sphere representation. 
As evident, here the average SOP are about the same for 
the different input polarization orientations, as expected. 




(a) (b) 



Fig. 6. The average polarization when varying the in- 
put wavelength for different input polarizations on (a) 
the original Poincare sphere and on (b) the modified 
Poincare sphere. 

5. Concluding remarks 

To conclude, we presented experimental and calculated 
results about the dynamics of modes and their SOP in 
multimode fibers. We show how the representation of 
these dynamics can be simplified so as to obtain a direct 
measure of the SOP at the output when knowing the 
input SOP. The results of this measure can then be used 
for determining the SOP for an arbitrary input SOP. 
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